Energy Shift and Wavefunction Overlap of Metal- Organic Interface-States 
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The properties of Shockley-type interface states between 7r-conjugated organic molecular layers 
and metal surfaces are investigated by time-resolved two-photon photoemission experiments and 
density functional theory. For perylene- and naphthalene-tetracarboxylic acid dianhydride (PTCDA 
and NTCDA) adsorbed on Ag(lll), a common mechanism of formation of the interface state from 
the partly occupied surface state of the bare Ag(lll) is revealed. The energy position is found to be 
strongly dependent on the distance of the molecular carbon rings from the metal and their surface 
density. Bending of the carboxyl groups enhances the molecular overlap of the interface state. 

PACS numbers: 73.20.-r, 78.47.-p, 79.60.Dp 
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The energetic position and wavefunction overlap of 
electronic states at the interface between layers of or- 
ganic molecules and metals is of fundamental interest 
for the design of organic semiconductor devices and for 
future applications of molecular electronics. Previous 
studies, both experimentally and theoretically, concen- 
trated on those electronic states that either result di- 
rectly from the chemical bonding at the interface or from 
the shift and broadening of localized molecular orbitals 
upon interaction with the metal substrateP" 8 -' However, 
also states intrinsic to metal surfaces are affected by ad- 
sorption and can become an important factor for the 
electronic coupling between metal and organic molecules. 
This has recently become apparent for 3,4,9, 10-perylene- 
tetracarboxylic-acid-dianhydride (PTCDA) on Ag(lll), 
a structurally very well characterized model system for 
the interface between 7r-conjugated organic molecular 
layers and metals™ A dispersing, free-electron-like elec- 
tronic state located 0.6 eV above the Fermi level was 
observed by tunnelling spect roscopy^ and by two-photon 
photoemission (2PPE) fi° * 11 llt originates from the Shock- 
ley state of the Ag(lll) surface which is upshiftcd from 
below the metallic Fermi level by as much as 0.7 eV. 
The initially partially occupied state becomes unoccu- 
pied and approaches the conduction band of the organic 
semiconductor.EOffil] 

In this Letter, we address the question which prop- 
erties of the molecule-surface interaction determine the 
energetic position of such Shockley-type interface states 
(IS) and which factors facilitate a large overlap of the 
state with the molecular layer in order to tailor the de- 
gree of electronic coupling. For rare-gas adsorbates on 
noble metals it is well established that their adsorption 
also leads to an upshift of the Shockley state, albeit by 
much smaller values, and that the shift increases system- 
atically from the weakly bound Ar to the more strongly 



interacting Xe! 14 * 15 ! The comparison of our 2PPE experi- 
ments for 1,4,5,8-naphthalene-tetracarboxylic acid dian- 
hydride (NTCDA) on Ag(lll) with PTCDA/ Ag(lll) 
shows an analogous trend. The structurally similar, but 
smaller NTCDA molecules cause a smaller IS upshift 
than the larger PTCDA molecules. The time-resolved 
2PPE data indicate a similar overlap of the wavefunction 
with the metal in both cases. Density functional calcula- 
tions explain the smaller shift by a larger distance of the 
NTCDA carbon plane from the surface. Surprisingly, 
in the metastable disordered PTCDA phase, when the 
molecules interact more strongly with the substrate™ we 
observe a smaller and not a larger IS upshift. Lifetime 
measurements indicate a reduced overlap of the wave- 
function with the metal. Theory relates this finding to 
the bending of the carboxyl groups which enables the IS 
wavefunction to extend more into the molecule. 

The experiments have been conducted under ultra- 
high vacuum conditions using a 2PPE setup described 
previously.^ The photon energies of the blue and ultra- 
violet (UV) laser-pulses were 3.10 eV and 4.70 eV with 
pulse lengths of 47 fs and 70 fs, respectively. The overall 
energy resolution of the setup was 65 meV, the angular 
resolution 1.2°. Thin organic films of temperature gra- 
dient sublimated PTCDA and NTCDA were deposited 
onto sputter/annealed, clean Ag(lll) surfaces at a rate 
of 0.2-0.4 ML/min. Thermal desorption of multilayer 
films resulted in a commensurate monolayer for PTCDA 4 
and the relaxed monolayer for NTCDAPS Sample cooling 
to 90 K was performed rapid enough (> 15 K/min) to 
suppress the inverse melting phenomenon for NTCDA. 18 
The PTCDA-LT phase was prepared by adsorption at a 
sample temperature of 90 K. 

The DFT calculations have been performed using the 
pseudopotential method as implemented in OPENMX 
code within the local density approximation for the 
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FIG. 1: Energy-, angle- and time-resolved 2PPE data of monolayer films of NTCDA, PTCDA on Ag(lll). (a) Normal emission 
(UV+UV) 2PPE spectra of commensurate PTCDA and NTCDA monolayers recorded at 300 K and 90 K as indicated, and of 
the disordered LT-phase of PTCDA at 90 K. The insets sketch the molecular structure for PTCDA (top) and NTCDA (bottom), 
(b) Angular dependence of the peak maxima of the interface state of the ordered PTCDA (squares, red) and NTCDA (diamonds, 
green) phases measured at 300 K and the PTCDA low-temperature phase (dots, blue) measured at 90 K together with projected 
Ag(lll) bulk bands (gray shaded area) and the Shockley state of the clean surface (black), (c) Excitation schemes employed to 
2PPE spectroscopy (left) and time-resolved 2PPE (right), (d) Time-resolved 2PPE data recorded at the peak maxima of the 
interface states of the different PTCDA and NTCDA monolayers. Lifetimes were extracted from best fits of a rate equation 
model using Gaussian laser pulses that reproduce the measured (Blue+UV)-cross-correlation on the clean surface (dashed line). 



exchange-correlation potential! 12 * 19 ! A slab containing a 
9-layer silver film was employed. The relaxed NTCDA 
monolayer was attached on one side of the film, the vac- 
uum region was chosen to correspond to 6 silver inter- 
layer spacings. The resulting periodic supercell has a 
rectangular unit cell of 11.56 x 15.02 A 2 and contains 
216 silver atoms and 48 atoms belonging to two NTCDA 
moleculesPSl The surface Brillouin zone (SBZ) was sam- 
pled with a 3 x 3 k-point mesh. 

Figs. l(a,b) display the results of angle-resolved 2PPE 
spectroscopy performed with a single-color (UV+UV)- 
excitation scheme (cf. Fig. 1(c)). In agreement with 
previous results, the well-ordered herringbone phase of 
PTCDA/ Ag(lll) exhibits a strongly dispersingstate 
~ 0.6 - 1.0 eV above the Fermi level (E F )^ For 
one monolayer the band minimum at k\\ = is 
E — E F = 0.57 ± 0.02 eV at room temperature. Its en- 
ergy increases slightly when additional layers are ad- 
sorbed or when the sample is cooled to liquid nitrogen 
temperatureP Similar to PTCDA, the smaller NTCDA 
molecules can be adsorbed in a commensurate monolayer 
on Ag(lll)P2l We observe an unoccupied dispersing in- 
terface state also for this system (green lines and data 
points in Figs. l(a,b)). It is located 0.2 eV closer to the 
Fermi level (E - E F = 0.38 ± 0.02 eV at jfcy = 0) than 
the IS of PTCDA. The effective masses at the T-point 
of both IS are similar (PTCDA: 0.46 ± 0.1 m e , NTCDA: 
0.47 ± 0.1 m e ) and both values are close to that of the 
Shockley state of the clean metaf^ (0.42 m e ). 



The results of the time-resolved 2PPE experiments 
(Fig. 1(d)) underline the similarities between the IS of the 
commensurate PTCDA and NTCDA layers. For these 
experiments we made use of a two-color (blue+UV)- 
excitation scheme. The unoccupied IS is populated by 
an intense blue pump pulse and the transient popula- 
tion is photoemitted by a weaker UV pulse with variable 
time delay (Fig. 1(c)) P The lifetime of 53 ± 3 fs deter- 
mined for the IS of the PTCDA herringbone phase at 
300 K and feii = agrees well with our previous results.^ 
For NTCDA we measure a considerably longer lifetime 
of 115 ±10 fs. 

Before comparing these two values we note that 
both IS are located in the gap of the projected bulk 
bands of Ag(lll). Electrons excited into these states 
will decay inelastically by electron-hole-pair creation 
at a rate determined by the wavefunction overlap 
with the metal, the screening of the Coulomb inter- 
action and the available phase space in the metal.^ 
Under the realistic assumption of a similar screened 
Coulomb interaction at both interfaces (NTCDA/Ag 
and PTCDA/Ag) and taking into account the different 
phase-space factors for the decay at the different energies 
above the Fermi level, the measured lifetimes of the 
IS allow us to compare their overlap with the metal. 
In fact, the ratio of the lifetimes tntcda/tptcda 
[116/53 = 2.19] is very close to the inverse phase 
space factor (-Eptcda - E F ) 2 /(E NTCDA - E F ) 2 
[(0.57) 2 /(0.38) 2 = 2.25] of a free electron gas close 
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FIG. 2: (a) Band structure of the NTCDA monolayer on the nine-layer Ag(lll) film in the V — > Xi direction of the surface 
Brillouin zone for the bend molecular geometry. Filled circles indicate wave vectors k, at which the spatial distribution of the 
IS wave function was analyzed, the thick line represents the parabolic approximation of the corresponding dispersion, arrows 
indicate IS energies at F for different geometries (from top to bottom: red: d = 2.85 A-fiat, orange: d = 2.997 A-flat, blue: less 
dense packing), (b) Spatial distribution of the IS probability density p xy at the F-point in the vicinity of the molecular plane 
for the flat (top) and the bent-down geometry (bottom), (c) IS probability density p z perpendicular to the surface, rry-averaged 
over the unit cell (top: p z for the three considered NFCDA geometries; bottom: difference Ap z between the bent and the flat 
geometry at the distance of 2.997 A. 



to the Fermi level. 23 This suggests a similar overlap 
of both states with the metal. We arrive at the same 
conclusion when we compare the lifetimes of the interface 
states of NTCDA and PTCDA along their dispersing 
bandsP At k\\ = 0.22 A" 1 , e.g., the NTCDA-IS has 
a similar energy as the PTCDA-IS at k\\ = 0.19 A -1 
(E - E F ~ 0.8 eV) and its lifetime (38 fs, Fig. 1(d)) is 
close to that of the PTCDA-I^l (33 f s )). 

Our results thus indicate that the commensurate 
NTCDA and PTCDA layers have a qualitatively very 
similar influence on the Shockley surface state of 
Ag(lll), except for a larger upshift induced by the more 
strongly interacting PTCDA. 

This picture changes when we turn to the results ob- 
tained for the low temperature phase of PTCDA. Ad- 
sorption at sample temperatures below 150 K leads to 
an incommensurate PTCDA monolayer in which the 
molecules form dendritic islands. The reduced lateral in- 
teraction of the molecules in this phase goes along with 
a stronger binding to the surfaced Surprisingly, this 
stronger binding is not accompanied by a larger IS up- 
shift as compared to the commensurate PTCDA phase. 
It is observed at an energy of 0.50±0.03 eVJ^i.e. 0.07 eV 
lower than the IS of the more weakly bound commensu- 
rate phase (Fig. 1(a)). The electron lifetime is 54 ± 3 fs 
(Fig. 1(d)). Since the IS lifetime of the commensu- 
rate PTCDA phase is strongly temperature dependent, 17 
we have to compare this value, obtained at 90 K, with 
lifetime measurements of the commensurate phases per- 



formed at the same temperature (Fig. 1(d)). The life- 
time of 54 ± 3 fs is almost twice the 26 fs measured for 
commensurate PTCDA at 90 K and even longer than the 
IS lifetime of NTCDA (43 fs at 90 K) though the latter 
state is closer to the Fermi level. Following the argu- 
ments above, the time-resolved experiments thus indicate 
that the wave function of the IS of the low-temperature 
PTCDA phase has less overlap with the metal than that 
of the ordered PTCDA and NTCDA layers. 

In order to obtain a microscopic understanding of the 
origin of the observed energy shifts of the IS and its 
wavefunction overlap, we performed DFT calculations 
for three different geometries of the NTCDA monolayer. 
First, the molecules stay in their flat gas phase geome- 
try, while they arc artificially kept at the experimentally 
determined adsorption distance of the carbon plane (d = 
2.997 A). Second, the flat molecules are moved slightly 
closer to the surface (d = 2.85 A), i.e. to the adsorption 
distance of the ordered PTCDA monolayer.— Third, the 
experimental adsorption geometry of NTCDA with bent 
carboxylic oxygen atoms is simulated (d — 2.997 A)P3 
Comparing calculations for one molecule (NTCDA) in 
different geometries allows us to pinpoint the relevant 
factors that determine properties of the interface state 
very clearly, although we can of course not expect full 
quantitative agreement with our experimental results of 
NTCDA and PTCDA molecules in different adsorption 
geometries. 

Fig. 2(a) shows the calculated band structure of 



NTCDA on Ag(lll) in the experimentally observed ge- 
ometry. Due to the presence of the molecular monolayer, 
the period of the unit cell has grown in comparison with 
the clean Ag(lll) surface, and the corresponding surface 
Brillouin zone has reduced its size. As a result, bands in 
the SBZ of the clean-surface unit cell get "folded" into 
the reduced SBZ. This leads to a surface band struc- 
ture that no longer exhibits the projected bandgap at 
the f -point (cf. Fig. 1(b)). As the (metal) electrons will 
actually have very small scattering probability at the bor- 
der of this SBZ, this fact does not change the arguments 
given above for the decay of electrons excited into the 
IS, but it complicates the search for the IS in the calcu- 
lated band structure. For that purpose we have analyzed 
the spatial probability distribution of states at different 
k-points along symmetrical directions, averaged over the 
spatial coordinates xy within the supercell (Fig. 2). 

Characteristic for the IS is a probability density along 
the surface normal p z that resembles that of the Shock- 
ley state of bare Ag(lll) in the silver film. It has a 
pronounced maximum in the region of the metal-organic 
interface and additional peaks around the center of the 
NTCDA molecules (Fig. 2(c)). The dispersion of the IS 
reflects its metallic character and can be described by a 
free-electron like parabola (Fig. 2(a)). Within our accu- 
racy, the effective mass is independent of the geometry 
(0.44 m e for flat and 0.45 m e for bent molecules) and close 
to that of the Shockley statePl 

In contrast to the effective mass, the IS energy position 
at the F-point is very sensitive to the exact geometry of 
the NTCDA adlayer (arrows in Fig. 2(a)). For the dis- 
tance of 2.997 A, the bending only leads to a slight upshift 
of 41 meV, from 243 meV to 284 meV. The reduction of 
the distance by merely 5% to 2.85 A, however, induces 
an upshift of 112 meV to 355 meV. The, as compared to 
NTCDA, higher IS energy of the commensurate PTCDA 
layer can therefore to a large extent be explained by its 
slightly smaller adsorption distance. 

Another important factor is the relative density of 
molecules per unit cell area in general and carbon rings 
in particular. Keeping only one NTCDA molecule in the 
unit cell, brings the IS energy ~ 100 meV closer to the 
Fermi level. This finding explains why, in spite of the 
shorter adsorption distance, the IS is observed at lower 
energy in the disordered, and consequently less densely 



packed LT-PTCDA phase. Furthermore, the previously 
calculated IS energy for PTCDA, with an adsorption dis- 
tance of 2.85 A, is 469 meV at the f-pointP Although 
this value, corresponding to the bent geometry, should 
not be quantitatively compared with the 355 eV, ob- 
tained for flat NTCDA at 2.85 A, the bending of the car- 
boxyl groups is unlikely to be the dominating factor for 
the larger upshift of the PTCDA-IS. Rather this compar- 
ison suggests that, together with the reduction of the dis- 
tance from the metal surface, the higher density of carbon 
rings per unit cell area can account for the experimen- 
tally observed energy difference of 200 meV between the 
commensurate NTCDA and PTCDA layers. Both factors 
steepen the potential at the metal/organic as compared 
to the metal/vacuum interface and thus cause the Shock- 
ley state of the metal to shift to higher energies. 

Whereas the bending of the carboxyl groups does not 
have a large effect on the IS energy position, it signifi- 
cantly changes the lateral probability distribution p xy of 
the wave function. In particular, the amplitude at the 
bridging oxygen atoms and at the carbon atoms of the 
carboxyl groups is greatly enhanced as compared to the 
flat molecule (Fig. 2(b)). In addition, p xy is redistributed 
within the naphthalene core. Averaged over the whole 
SBZ, p z becomes larger on the molecule and diminished 
in the metal (Fig. 2(c)). This explains the longer life- 
time of the IS of LT-PTCDA compared to the ordered 
phase since NIXSW experiments have shown a signifi- 
cantly stronger bending of the carboxyl groups of these 
molecules towards the metal™ 

In conclusion, we investigated the formation mecha- 
nism of unoccupied Shockley-derived metal-organic inter- 
face states. With adsorption of the PTCDA and NTCDA 
molecules on Ag(lll) the Shockley state of the metal is 
shifted to energies above the Fermi level. The new en- 
ergy position is primarily influenced by the presence of 
a 7r-electron system. It depends on adsorption distance 
and areal concentration of carbon rings. In contrast, the 
spatial probability density of the state, which has a pro- 
nounced maximum between the terminating silver atoms 
and the carbon plane, is mainly determined by the molec- 
ular geometry. Bending of the carboxylic end groups re- 
duces the metallic and enhances the molecular overlap. 
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